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Abstract — Spatially-coupled low-density parity-check codes at- 
tract much attention due to their capacity-achieving performance 
and a memory-efficient sliding-window decoding algorithm. On 
the other hand, the encoder needs to solve large linear equations 
to terminate the encoding process. In this paper, we propose 
modified spatially-coupled codes. The modified (di, d T ,L) codes 
have less rate-loss, i.e., higher coding rate, and have the same 
threshold as {d\,d T ,L) codes and are efficiently terminable by 
using an accumulator. 

I. Introduction 

Spatially-coupled (SC) low-density parity-check (LDPC) 
codes attract much attention due to their capacity-achieving 
performance and a memory-efficient sliding-window decoding 
algorithm. The studies on SC-LDPC codes date back to the in- 
vention of convolutional LDPC codes by Felstrom and Zigan- 
girov [1]. They introduced a construction method of (d\,d T )- 
regular convolutional LDPC codes from (d\, <i r )-regular block 
LDPC codes. The convolutional LDPC codes exhibited better 
decoding performance than the underlying block LDPC codes 
under a fair comparison with respect to the code length. 
Lentmaier et al. observed that (4,8)-regular convolutional 
LDPC codes exhibited the decoding performance surpassing 
the belief propagation (BP) threshold of (4,8)-regular block 
LDPC codes [2]. Further, the BP threshold coincides with the 
maximum a posterior (MAP) threshold of the underlying block 
LDPC codes with a lot of accuracy. Constructing convolutional 
LDPC codes from a block LDPC code improves the BP 
threshold up to the MAP threshold of the underlying codes. 

Kudekar et al. named this phenomenon "threshold sat- 
uration" and proved rigorously for the binary-input erasure 
channel (BEC) [3] and the binary-input memoryless output- 
symmetric (BMS) channels. [4]. In the limit of large d\,d r ,L 
and w, the SC-LDPC code ensemble (d\, d T , L, w) [3] was 
shown to universally achieve the Shannon limit of BMS 
channels under BP decoding. This means the transmitter does 
not need detail statistics of the channel but needs to know 
only the channel capacity. Such universality is not supported 
by other efficiently-decodable capacity-achieving codes, e.g., 
polar codes and irregular LDPC codes. According to the 
channel, polar codes need frozen bit selection and irregular 
LDPC codes need optimization of degree distributions. We 
note that recently Aref and Urbanke proposed SC rateless 
codes [5] which are conjectured to universally achieve the 
capacity of BMS channels without knowing even the capacity 
of the channel at the transmitter. 

SC-LDPC codes are constructed from L LDPC block codes 
and have lower coding rate than the underlying block codes. 



As increasing coupling number L, the rate-loss could be 
decreased and the block error rate are degraded. Kudekar et 
al. proposed to mitigate the rate-loss in [6]. 

Pusane et al. proposed an efficient encoding and decoding 
method of convolutional LDPC codes in [7]. Encoding of 
convolutional LDPC codes can be divided into two stages: 
sequential encoding and termination. The sequential encoding 
process can be efficiently calculated with computational cost 
O(LM), where M is the lifting number for protograph codes. 
On the other hand, termination process needs a multiplication 
of dense matrix of size, which needs computational cost 
0(M 2 ). Since M and L are chosen so that L <C M, termi- 
nation needs more computations than decoding and sequential 
encoding. 

In this paper, we propose modified spatially-coupled codes. 
The modified (d\,d I: L) codes have less rate-loss, i.e., higher 
coding rate, and have the same threshold as (d\,d T ,L) codes 
and are efficiently terminable by using an accumulator. 

II. Preliminaries 

In this section, we briefly review (d\,d r ,L) codes intro- 
duced by Kudekar et al. [3]. We assume 4 =: i; £ Z and 
k > 2. 

The SC-LDPC codes are defined by the following proto- 
graph codes [8]. The adjacency matrix of the protograph is 
referred to as a base matrix. The base matrix of (d\, d T , L) code 
is given as follow. Let H(d\, d T , L) be an (L + d\ — 1) x kL 
band binary matrix of band size d r x d\ and column weight d\, 
where the band size is height x width of the band. We refer 
to L as coupling number. For example, H(4, 12, 9) is given 
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#(4,12,9) 



The protograph of the base matrix H (4, 12, 9) is given in 
Fig. |3] The protograph of (d\, d r , L) codes have kL variable 
nodes and L + d\ — 1 check nodes. The degree of variable 
nodes are all of degree d\ = 4. On the other hand, the degree 
of check nodes are not uniform, check nodes at center are of 
degree d T = 12 and check nodes near boundaries have lower 
degree. 
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Fig. 1. Parity-check matrix H(d\ = 3, d T = 6, L = 9, M) of (di,d T ,L) codes. Each Pj j represents an M X M random permutation matrix. The red 
submatrix is used for termination. 
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Fig. 2. Parity-check matrix H(d\ = 3, d T = 6, L = 9, M) of modified (di, d r ,L) codes. The red submatrix is used for termination. 



The Tanner graph of (d\, d T , L) code is obtained by making 
M copies of protographs of H(d\, d Yl L) and connecting edges 
among the same edge types. The parameter M is referred to 
as lifting number. The parity check matrix H(d\, d T , L, M) of 
a (di, d T , L) code is given by replacing each 1 in H {d\, d T , L) 
with an M x M random permutation matrix and each with 
an M x M zero matrix. Let Pij for i £ [1,L + d\ — 1], 
j G [1, dr] be a binary random M x M permutation matrices. 
An example H(4, 12, 9, M) is given in Fig. [TJ 

Denote a code word of (di,d T ,L) code by [x x , . . . , x kL ), 
where x ■ e F^ 7 for j G [1 , fcL] . Each M rows of the i-th 
parity-check equations for H(d\, d T , L, M) is written as 

,/.-, II (1) 

for i = 1, . . . , L + di — 1. We assumed Xj = for j ^ [1, /cL] 
for simplicity of notation. 

When the parity-check matrix H of a protograph code is 
full -rank, the coding rate is given by 1 — — , with n variable 
nodes and m check nodes. Unless otherwise specified, we 
assume parity-check matrices are full-rank. Hence, the design 
coding rate of (d\,d r ,L) codes is given by 

R = kL-(L + dj- 1) = k-l_ _ dj-1 
kL k kL 

The rate-loss from the coding rate (k — l)/fc of (d\, d r ) codes 
is ^r-i which vanishes as 0(1/ L). 

Lentmaier et al. [9] observed that as L increases, the 
BP threshold values e BP of (d\,d Ti L) codes approach the 
MAP threshold value of (d\, d T ) LDPC codes. Large coupling 
number L is not preferred since the block error rate of SC- 
LDPC codes tends to be degraded when L is large. Hence, 



mitigating the rate-loss while keeping the BP threshold is 
desired. 

III. Issues in Encoding of Spatially-Coupled Codes 

Encoding convolutional LDPC codes has been investigated 
by Pusane et al. in [10]. Encoding of (d\, d r , L) codes involves 
two parts, i.e., sequential encoding and termination. 

For simplicity, we assume M = 1 unless otherwise speci- 
fied. Since there are kL variable nodes and L + d\ — 1 check 
nodes, the number of information bits iVi n f is given by 

N inio = kL - (L + di-l). 

In the sequential encoding stage, for each section i = 1, . . . , 
the parity bit x ik is determined by the z-th check-node 
constraint ([TJ and by using coded bits x 1} . . . ,X/i_i\j c which 
was determined in the previous stages and k — 1 information 
bits • • • ,%_ik-i- The total number A^ seq of parity 

bits which are determined in the sequential encoding stages is 
given by 
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For the example of H(A, 12, 9, M), N soq (4, 12) = 7. We need 
to determine the remaining N ieTm undetermined parity bits 
from the remaining A^ n fo— (k— l)N soq information bits, where 

iVterm = kL - {N info + N scq ) = d\ - 1 + ^— j . 

For the example of H(4, 12, 9, M), N teim (4, 12) = 5. 

In Fig. [3] protograph of (4, 12, 9) code is shown. Gray nodes 
represent information bits. Light red nodes and red nodes are 
parity nodes determined in the sequential encoding stage and 
termination stage, respectively. 



A. Sequential Encoding Stage 

Let us see how the i-th parity bit is determined in the 
sequential encoding stage. Define syndrome at section i as 
follows. 



For the example of H (4, 12, 9, M), we need to solve linear 
equations H teim (x 23 , ■ • ■ ^ 27 ) T = (s 8 , . . . , s 12 ) T with 
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where, for simplicity of notation, we define Xj = for 
j [1 , kL] . We start from s x = 0. In sequential encoding 
part, for each section i = 1,2, . . ., one parity bit node x ki 
is sequentially determined from k — 1 information bit nodes 
■ ■ • ,x ki _ 1 and syndrome s { . From (U} and (O, we 

have 

2i+ E PiMik-dr+j = Q- W 
From this, x fei can be determined as follows. 



PijZLik-dt+j ) 



(5) 



d r -l 

E 

j =c ; r -fc+i 

Since P^ 1 is an M x M permutation matrix, (0 can be 
accomplished with O(M) computational costs. In Fig. [3] 
protograph of (4, 12, 9) code is shown. Sequential encoding is 
accomplished sequentially from left for each fc(= 3) variable 
nodes. First, x 3 is determined from Next, x 6 is deter- 



mined from s 



21 2L4 



ic. Then x Q is determined from s 



31 ±7: ±8' 



This continues until x 21 is determined from s 7 ,x lg ,x 2 o- The 
sequential encoding process continues until the number of the 
remaining information bit nodes reaches less than k — 1. After 
that, we move to termination stage. 

B. Termination 

In the termination stage, we determine the remaining iVt erm 
parity bits so that parity-check equations ([T) are satisfied for 
i E [N scq + 1 , -ZV soq + A^term] ■ We define the i-th syndrome 

Sj (i G [A soq + 1, N se q + iVterm]) by the contribution from The rate-loss from the rate 
determined bits to the i-the check nodes. To be precise, 

E 

]':=ik-d T +j<Ni„ lo +N s , 

dr-ik+Ninfv+Nscq 



(7) 



In Fig. [B ffterm is written in red in H(4, 12, 9, M). In 
FigOred nodes represent bit node involved in the termination 
process. 

IV. Modified Spatially-Coupled Codes 

In this section, we propose modified (di,d T ,L) codes. We 
give an explanation how the modified codes can be efficiently 
terminated. 

A. Modified (d\,d T , L) Codes 

The base matrix H(d\,d T , L) of modified (d\,d r ,L) codes 
is obtained by removing d\ — 2 bottom rows of the base matrix 
H(d\,d r , L) of (di,d T ,L) codes. For example, H(A, 12,9) is 
given as follows. 

rill 
lllili 
111111111 
111111111111 
111111111111 
111111111111 
111111111111 
111111111111 
111111111111 
llllillll. 
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The protograph of modified codes is obtained by removing 
(d\ — 2) right most check nodes and their incident edges from 
the protograph of (d\, d T , L) codes. In Fig. |H for example, we 
show protograph of modified (d\ = 4, d r = 12, L = 9) codes. 
The coding rate of the modified (d\, d T , L) codes is give by 



R(d h d T ,L) 



kL -(L + l) = k-1 _ J_ (g) 
kL k kL 

of (di,d T ) codes is 
Comparing (O and dHJ, we can see that rate-loss is mitigated 
as (d\ — l)-times as much by modifying. 
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Pi,j-Lik-d T +j- 



V. Efficient Encoding of Modified Codes 

In this section, we explain how modified (d},d T ,L) codes 
can be efficiently encoded. 



Denote the right bottom N teTm M x N te rmM submatrix 
of H(di,d T ,L,M) by i?term- The termination process of 
(d\,d T ,L) codes is equivalent to solve the following equation 



Hi 



iS-L+dl-Nu, 



+D • ■ • 7 SLicl) 
m> • • - ,s L+dl _ 



A. Sequential Encoding 

Sequential encoding is performed in the same way as 
(d\, d T , L) codes. The number N- ln f of information bits of 
modified (d\,dr,L) codes is given by 



(6) 



fo 



kL-(L+l). 



When -ffterm is not full-rank, one can modify some entries 
of i/term so that H torm becomes full-rank and the decoding 
performance remains almost the same. The inverse of i/term is 
not sparse in general. Hence, one needs 0(M 2 ) computational 
cost to solve the linear equations ((6). 



The total number iV se q of parity-bits which are determined in 
the sequential encoding stages is given by 



N seq 



info 



l - 1 



2:1 £.2 —3 —4 —5 —6 —7 —8 —9 —10—11—12 —13—14—15—16 —17—18—19—20—21—22^23—24—25—26—27 

Fig. 3. Protograph of (dj = 4, d r = 12, L = 9) codes. Gray nodes represent information bits. Light red nodes represent parity bits. Red nodes represent 
parity bits calculated in termination. Sequential encoding is accomplished sequentially from left for each k(= 3) nodes. Termination requires to solve linear 
equations involving as much as five variable nodes x 23 , . . . ,2:27 snown m 0- One needs 0(M 2 ) computational cost to solve the linear equations. 




Fig. 4. Protograph of modified (d\ = 4, d T = 12, L = 9) codes. This protograph is produced by removing (d\ — 2) right most check nodes from the 
protograph of (d\ = 4, d T = 12, L = 9) codes. Termination is accomplished by solving linear equations involving only two variable nodes 
One can solve the linear equations by an accumulator with 0(M) computational costs like RA codes. 

TABLE I 

Comparison of coding rate and BP threshold values for modified {d\, d r ,L) codes and (di, d T , L) codes. 



(d h dr,L) 


? BP 
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(dl,d T ,L) 


-BP 


f BP 
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R 


(3,6,9) 


0.49174 


0.51203 


0.44444 


0.38889 


(3,9,9) 


0.32157 


0.33305 


0.62963 


0.59259 


(3,6,17) 


0.48816 


0.48876 


0.47059 


0.41177 


(3,9,17) 


0.31997 


0.31995 


0.64706 


0.62745 


(3,6,33) 


0.48815 


0.48815 


0.48485 


0.46970 


(3,9,33) 


0.31965 


0.31965 


0.65657 


0.64647 


(3,6,65) 


0.48815 


0.48815 


0.49231 


0.48462 


(3,9,65) 


0.31965 


0.31965 


0.66154 


0.65641 


(4,8,9) 


0.50158 


0.51938 


0.44444 


0.33333 


(4, 12,9) 


0.33282 


0.33282 


0.62963 


0.52963 


(4,8,17) 


0.49774 


0.49787 


0.47059 


0.41177 


(4, 12, 17) 


0.33025 


0.33033 


0.64706 


0.60784 


(4, 8, 33) 


0.49774 


0.49774 


0.48485 


0.45455 


(4,12,33) 


0.33025 


0.33025 


0.65657 


0.63636 


(4,8,65) 


0.49774 


0.49774 


0.49231 


0.47692 


(4, 12,65) 


0.33025 


0.33025 


0.66154 


0.65128 



B. Termination 

It follows that the number of parity bit nodes we need to 
determine in termination process is 



= kL- N info - N scq = 2 



Note that iVterm( = 2) does not depend on d\ and d r . This 
means that we always need deal with only two parity bits 
to terminate the encodin process. In the termination stage, we 
determine the remaining Nterm(= 2) parity bits so that parity- 
check equations (H) are satisfied for i 6 {TVseq+1, A r soq +2} = 
{L,L + 1}. 

We define the i-th syndrome s, l e F^ 1 (i 6 {L,L+ 1}) by 
the contribution from determined bits to the i-the check nodes. 
To be precise, 



dr-ik+N in t +N e , 

E 



Pi. 



i,j—ik—d T +j 



-ik+kL-2 

E - 



i,j3Lik-d r +j- 



Denote the right bottom iVterm-^ x N ieim M submatrix of 
H(d\,d r , L, M) by -ffterm- The termination process of modi- 
fied (d\,d T ,L) codes is equivalent to solve the following linear 
equation. 

Htorm{x kL _ 11 x kL ) T — (S L ,S L+1 ) T , (9) 

Ph,d,k 



-l:±kLJ 

PL,d } k-l 
Pt. 



L+l,(d!-l)fc-l PL+l,{di-l)k 

In Fig. 12 the submatrix Htcrm is printed in red. In Fig. |4] 
parity bits involved in termination are drawn in red 

Since the sum of all rows of -fftcrm is € F| Ai , H x 
is not full-rank. We need a slight modification to ensure that 
S-kL-i an d Z-kL can t> e determined from s L and S L+1 . We 
modify H tC i- m as follows. 

Im I'm 
Iai Im 



im 

2 ' 11 term 



(10) 




Fig. 5. Comparison of convergence speed of decoding erasure rate pf and pf 1 at the £-lh BP iteration at section i S [1, 27] for (4,12,9) codes (left) and 
modified (4,12,9) codes (right) BEC(e=0.3) 
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where Im is the identity matrix of size M x M. The effect for 
decoding performance from this modification can be negligibly 
small when lifting number M is large. The equation (O with 
-Htcrm ( flOb can be solved by 2M accumulations with an 
accumulator in the following way. Let us define x kL _ 1 =: 
(xkL-i,i, ■ • ■ , XkL-i,At), x kL =; (xkL,i, ■ ■ ■ , XkL-i,M)- The 
solution x kL _ 1 and x kL1 in the equation © can be solved by 
sequentially calculating 

XkL-l.i = XkL.i-1 + s L,i 
XkL,i = XkL-l,i + s L+l,i 

for i = 1 . . . , M, where we assumed x^l-i = 0. Thus, the 
modified (di,d T ,L) codes can be terminated with computa- 
tional costs 0(M). 

VI. BP Threshold 

In this section, we give comparison of BP threshold values 
for (di,d T ,L) codes and modified (d\,d T ,L) codes over the 
binary erasure channel (BEC) with erasure probability e. 

Table Q] shows BP threshold values e BP (c?i, d T , L) 
(resp. e BP (di, d T , L)) coding rate R(d\,d T ,L) 
(resp. R(d\, dr, L)) of (d\,d r ,L) codes (resp. modified 
(d\,d T ,L) codes). e sha represents the Shannon threshold for 
coding rate 1 — d\/d r . It can be seen that modified (d\, d T , L) 
codes have higher coding rate and the same threshold as 
(d\,d T ,L) codes within the 5 digits precision. 

Let p\ and p\ be the decoding erasure rate of x± at the 
£-th BP iteration for (4,12,9) codes and modified (4,12,9) 
codes, respectively. Figure [5] shows how pf' and pf con- 
verge to 0. It can be seen that modified codes have slower 
convergence. 



VII. Conclusion 

We propose modified spatially-coupled codes. The modified 
codes have less rate-loss, and have the same threshold as 
(d\,d r ,L) codes and are efficiently terminable by using an 
accumulator. 

Future works include analysis of convergense, finite length 
performance, and performance over other channels like [11] 
and [12]. 
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